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DEFORMATION AND INSTABILITY 
OF UNDERTHRUSTING LITHOSPHERIC PLATES 

ABSTRACT 

Recent confirmation of the relationship between seismology and 
tectonic plate movement has focused attention on the deformation 
and instability of lithospheric plates. In this paper we construct 
models of the underthrusting lithosphere for the calculation of dis- 
placement and deflection. The first of our three contributions is a 
mathematical theory that rigorously demonstrates the elastic in- 
stability in the decending lithosphere. The theory states that litho- 
spheric thrust beneath island arcs becomes unstable and suffers 
deflection as the compression increases. Thus, in the neighbor- 
hood of the edges where the lithospheric plate plunges into the as- 
thenosphere and mesosphere its shape will be contorted. The 
second contribution is our calculation of the lateral displacement. 

We show that, before contortion, the plate will thicken and contract 
at different positions with the variation in thickness following a 
parabolic profile. Finally, we explain the depth distribution of the inter- 
mediate and deep focus earthquakes in terms of plate buckling and 
contortion. 
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DEFORMATION AND INSTABILITY 


OF UNDERTHRUSTING LITHOSPHERIC PLATES 
INTRODUCTION 

The inclined zones of seismicity which characterize all active island-arc 
systems are curved along their length. For some cases they can be recognized 
down to about 600 km. Since most arcs exhibit a more or less well-defined zone 
of maximum and minimum in seismic activity in the mantle, it is reasonable to 
ask whether such features of seismic activities are associated with the deforma- 
tion and instability of the sinking lithosphere (Isacks, Oliver and Sykes, 1968; 
Roman, 1970). In this paper, we shall investigate the regions of instability and 
analyze the patterns of deformation of the underthrusting lithospheric plates. 

The principal difficulty we meet, when the problem of deformation of an unstable 
plate is to be posed in its mechanical form, consists in the behavior of the plate 
material. In order to obtain the first approximation, we assume that the plate 
is brittle enough to produce earthquakes by local fracture and elastic to main- 
tain its shape at a depth of 600 km after passing through perhaps 1000 km of mantle 
(Sykes, 1966; McKenzie, 1969). The other difficulty is of a mathematical nature 
and in order to find the solutions of the equations of the problem the simplest 
schemes and elementary methods of analysis are used. Solutions of the theory 
of plates are given in forms which facilitate the computation of bending, twisting, 
compression, contortion and displacement of the plate. For the present purpose, 
the thermal effect on the distortion of plate is neglected. We shall, also, discuss 
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the results of this analysis in comparison with the observed depth ranges (Fig. 
1) of seismic activities. 


STRAIN ENERGY 


We consider the plate element dxdy shown in Fig. 2 subjected to the bending 
moments M x , M y and twisting moments M xy and M yx . The strain energy due 
to these moments is given by (Love, 1944; Timoshenko and Woinowsky-Krieger, 


1959) 


dU, = -M x dy 


+ "M xy dy 


/ d 2 V 

y9x 2 

d 2 v 


dx) +(-M y dx) 


1 


9y9x 


dx I + I -M yx dx 


a 2 v 

^ y 2 

' d 2 V 
, 9x9y 


dyj 

dy) 


( 1 ) 


where in each term the first bracket gives the average total moment on the ele- 
ment and the second gives the angle through which it acts. 


The bending and twisting moments are defined by 


M x 

My 


Eh 3 / 9 2 r? 9 2 t;\ 

12(1- m 2 )\9x 2 M 9y 2 / 

Eh 3 / 9 2 r? 9 2 r?\ 

12(1 - M 2 ) Uy 2 M 9 x 2 / 


M 


xy 


Eh 3 

12(1 -m 2 ) 


(1 -m) 


9 2 t? 

9x9y 



( 2 ) 


where E is the modulus of elasticity, h is the thickness of the plate, p is the 


Poisson's ratio and r? is deflection. 
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By substituting equation (2) into equation (1) , we obtain the expression of 


this energy stored in the plate in terms of the deformation pattern. This is 


u ‘ 4 d // 


+ — 2 ( , - ii) 

| ^3x 2 3y 2 / 


3 2 t? 3 2 t? 
3x 2 3y 2 


v 3x3y . 


dxdy (3) 


where 


12(1 - m 2 ) 

Since the underthrusting litho spheric plate cannot sink through the denser mantle , 
itis subjected to abuoyancyT from the mesosphere. ThebuoyancyT is a compression 
load. Therefore the strain energy in the center plane of the plate due to deformation 


in the z direction is 

° 2 = ^ 77 ' t ( i ?) 2 dxdy (4) 

Adding equations (3) and (4) gives the total energy in the plate due to bending 
strain, twisting strain and strains in the center plane of the plate 



Equation (5) will be used for the analysis of plate contortion. 


CONTORTION 

The mode of plate contortion can be obtained by using the energy method. 
Let us consider a rectangular underthrusting lithospheric plate under bending, 
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twisting and compression. The three dimensions of the plate are a, b and has 
shown in Fig. 3a. A resultant component p of the forces due to thermal convec- 
tion in the mantle (Fig. 3b) is assumed to act at (£ , f ). Assuming simple support 
on the edges of the plate in the mesosphere and lithosphere makes it possible 
to choose a deflection in the form (Timoshenko and Woinowsky-Krieger, 1959) 

oo oo 


m7TX nwy 

2_y A mn Sin — Sm~ 
m = 1 n = 1 

Substitution of equation (6) into equation (5) yields 


U = 


7r 4 ab 

8 


m n 



( 6 ) 


ab 


m 

m n 


A 2 

"mn 


/ m7r\ 

V a / 


(7) 


The change in potential energy due to p is equal to 

W = p E 2 A mn Sin^ Sin^ 
m n a b 


( 8 ) 


The total energy in the system must then be minimized with respect to the 
amplitude term A mn , or 

8 (U- W) 


8A : 


0 


mn 


(9) 


By substituting equations (7) and (8) into equation (9), we have 


7r 4 ab 


DA mn [ + 


b 2 


-V + ^A mn T (SLY 


c . mTT% n7rf 
- p Sin Sm— — = 0 


( 10 ) 
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from which 


4pSin^Sm^ 
a b 
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7r 4 ab D 


/ m 2 

+ n 2 ^ 

I 2 + 1 


w 

+ b 2 "; 

7T 2 D 

a 2 


Putting equation (11) into equation (6), the deflection equation is 


( 11 ) 


V = 


4p 


Si^SinHi 


7r 4 abD m 


• S 2 — 

m n / m 2 + n 2 \ 

\ a 2 b 2 / 


+ 


rm 


7T 2 D 


T 


SinJEEsin^ 

a b 


( 12 ) 


In equation (12), T is a compression load, its sign will be negative and the 
value of this negative term may become great enough so that the denominator of 
equations (11) and (12) approaches zero. Under these circumstances, the de- 
flection 7 ? approaches a very large value irrespective of how small the 
normal force p may be. A possible mode of plate contortion for m = 2 and 
n = 1 is shown in Fig. 4. In Fig. 4, A and B are regions of maximum curvature 
and C is the neutral position. T will increase when the plate continues to thrust 
into the mesosphere. As T increases to a value such that the stresses in the 
plate approach their critical value, buckling may occur in region A or B. No 
fracture will appear in the neutral region C in which the contortion of the plate 
changes sign. Therefore, the simplest explanation, one direct consequence of 
the model of Fig. 4, is that the intermediate maxima in seismic activity for is- 
land arcs are restricted to the regions of instability A and B in the a stheno sphere. 
For the purpose of comparison, we have plotted the depth ranges of maxima in 
the seismic activity for several island arcs in Fig. 4. It shows clearly that the 
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depth ranges of maxima in seismic activity are correlated with the computed re- 
gions of instability. In a worldwide composite plot of seismic activity versus 
depth (Isacks, Oliver and Sykes, 1968), a minimum in activity appears near the 
depth of about 300 km. Thus, the existence of a neutral region in the underthrust- 
ing lithospheric plates seems to provide a geophysical explanation for the absence 
of seismic activity at this depth. 


LATERAL DISPLACEMENT 

In order to compute the lateral displacement, we consider a plate which is 
dragged against the mesosphere because of convection. To simplify the problem 
we assume that the deformation of the plate corresponds to a plane strain before 
contortion occurs. Fig. 5 represents an underthrusting lithospheric plate with 
thickness h and width b. The strain is considered in the xz plane. The compo- 
nents of displacement u and v of a point in the plate are the solutions of the follow- 
ing equations (Love, 1944; Timoshenko and Woinowsky-Krieger, 1959) 

+ — ) u + + iA) = o 

\0X 2 3z 2 / G \3x 2 3x3z/ 

/ii + iiV + A ± G/_9^U + = 0 

\3x 2 9z 2 / G \3x3z 3z 2 / 

where 

x ME 

(1 + ju)d - 2 M ) 

and G is the modulus of rigidity. In equation (13), body forces are neglected. 
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Differentiating these equations, the first with respect to x, the second with 


respect to z, and adding them together, we find 


(— + Ji) (is.* = o. 

\3x 2 3z 2 / \ 3x 3z/ 

3u 3v 

Therefore, the volume expansion + satisfies the Laplace equation. 


(14) 


The solutions of equations (13) are 

X + G 


u = - 


V = - 


2G 

X + G 
2G 


(— + — ^ 

\3x 3z / 

f 3u + 3v\ 
\3x 3z/ 


+ 0 , 


+ 0 , 


(15) 


where 0 t and 02 are harmonic, i.e. , 



The general solutions of equations (14) and (16) can be expressed in the 
following forms 

^ = A 0 + Z [A k (x + iz) k + B k (x - iz) k ] 

3x 3z u k = i L J 

0j = C 0 + Z [C k (x + iz) k + D k (x - iz) k ] (17) 

k = 1 

0 2 = E 0 + ^ [Er (x + iz) k + F k (x - iz) k ] 

where the coefficients A 0 , C 0 , E 0 , A k , B k , C k , D k , E k and F k are to be 
determined bythe stresses on the boundary of the plate. The stress -displacement 
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relationships are 


ff xx 


a 


zz 



( 18 ) 


By prescribing a xx , a xz and o zz on the plate boundary, we obtain 

v 2 

displacement v in equation (15) on the surface of the plate for “ « 


the vertical 
1 . 


v 


X 


” 2G(X + G) 




+ 


3X + 2G 
48X 



T 


09) 


Equation (19) shows that the plate, being thickened for small values of x, 
may form symmetrical buckles on both surfaces z = ±h/2. The displacement v 
vanishes for 


x = x„ = '/ 4 (b 2 + 3X h 7 


Vi 


( 20 ) 


For x > x 0 , v changes sign. That is, in this region, the plate begins to be 
contracted. This region corresponds to a phase change in the material (Anderson, 
1967). Therefore, the leading parts of the underthrusting' lithospheric plates may 
be deformed with a parabolic profile and the buckling in this unstable region seems 
to be the cause of the deep seismic activities (Fig. 5). It should be noted that, if 


we take into account the real case of boundary conditions and the effects of ther- 
mal stresses and body forces on the distortion of plate, there will be no 
symmetry with respect to the axis x or the buckling. However, it may be expected 
that the plate will show a tendency to contort as a whole and to thicken in the 
mesosphere. 


c 
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CONCLUSION 


We have developed mathematical models of the underthrusting lithospheric 
plates to determine the regions of instability and to compute the patterns of de- 
formation. Two types of plate buckling in the asthenosphere and mesosphere 
are discussed. The absence of seismological events at the depth of 300 km is 
explained in terms of the mode of plate contortion. In comparison with the 
worldwide composite plot of seismic activity versus depth, we conclude that the 
deep and intermediate - depth earthquakes are restricted to those regions of 
instability in which buckling of the plate may occur. 
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Figure 4. Mode of plate contortion and the depth ranges of intermediate maxima in seismic activities for 
several island arcs. The neutral region C of the plate corresponds to the depth of minimum in 
seismic activity. A and B are regions of instability in asthenosphere. 
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Figure 5. Plate buckling in mesosphere and the depth range of deep maxima in seismic activities for several 
island arcs. The leading part D of the plate is the region of instability near the boundary of 
mesosphere and asthenosphere. 


NASA-GSFC 



